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Abstract. We consider a simple scalar reaction-advection-diffusion equation with ignition- 
type nonlinearity and discuss the following question: What kinds of velocity profiles are 
capable of quenching any given flame, provided the velocity's amplitude is adequately large? 
Even for shear flows, the answer turns out to be surprisingly subtle. 

If the velocity profile changes in space so that it is nowhere identically constant, (or if 
it is identically constant only in a region of small measure) then the flow can quench any 
initial data. But if the velocity profile is identically constant in a sizable region, then the 
ensuing flow is incapable of quenching large enough flames, no matter how much larger is 
the amplitude of this velocity. The constancy region must be wider across than a couple of 
laminar propagating front- widths. 

The proof uses a linear PDE associated to the nonlinear problem and quenching follows 
when the PDE is hypoelliptic. The techniques used allow the derivation of new, nearly 
optimal bounds on the speed of traveling wave solutions. 



1. Introduction 

We consider a mixture of reactants interacting in a region that may have a rather compli- 
cated spatial structure but is thin across. A mathematical model that describes a chemical 
reaction in a fluid is a system of two equations for concentration n and temperature T of the 
form 

v 2 

(1) T t + u-VT = kAT + —g(T)n 

K 

2 

. Vn /rn\ 

n t + u ■ vn = — An g(T)n. 

Le k 

The equations (|l|) are coupled to the reactive Euler equations for the advection velocity 
u(x,y,t). Two assumptions are usually made to simplify the problem: the first is a constant 
density approximation || that allows to decouple the Euler equations from the system ([!]) and 
to consider u(x,y,t) as a prescribed quantity that does not depend on T and n. The second 
assumption is that Le = 1 (equal thermal and material diffusivities) . These two assumptions 
reduce the above system to a single scalar equation for the temperature T. We assume in 
addition that the advecting flow is unidirectional. Then the system ([]]) becomes 

(2) T t + Au{y)T x = kAT + —f(T) 

K 

T(0,x,y) = T (x,y) 

with f(T) = g(T)(l — T). We are interested in strong advection, and accordingly have written 
the velocity as a product of the amplitude A and the profile u(y). In this paper we consider 
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nonlinearity of the ignition type 

(i) f(T) is Lipschitz continuous on < T < 1, 
(3) (ii) /(l) = 0, 36 such that /(T) = for < T < O , f(T) > for T > O , 

(hi) f(T)<T. 

The last condition in ([3]) is just a normalization. We consider the reaction-diffusion equation 
(fj) in a strip D = {x G R, f/ £ [0, H]}. Equation (§) may be considered as a simple model 
of flame propagation in a fluid 0, advected by a shear (unidirectional) flow. The physical 
literature on the subject is vast, and we refer to the recent review |f22| for an extensive 
bibliography. The main physical effect of advection for front-like solutions is the speed-up of 
the flame propagation due to the large scale distortion of the front. The role of the advection 
term in (0) for the front-like initial data was also a subject of intensive mathematical scrutiny 
recently. Existence of unique front-like traveling waves has been established in ||, and 
their stability has been studied in [IS], [19], [21]. A traveling front is a solution of (0) of the 



form 

(4) T(t,x,y) = U(x-c A t,y) 
such that 

(5) lim U(s,y) = l, lim U(s,y) = 0, U s (s,y)<0. 



-oo s— >+oo 



The monotonicity property is not required for traveling wave solutions, but it is always present 
in the situation we consider. The speed-up of the fronts by advection mentioned above may be 
quantified as the dependence of the traveling front speed ca on the amplitude A. Variational 
formulas for ca were derived in [[II], |TT| . The latter work also contains results on the asymptotic 
behavior of ca when A is small for some classes of shear flows, as well as upper bounds on ca 
linear in A. An alternative approach to quantifying advection effects was introduced by the 
present authors in |7| [LSI . It is based on the notion of the bulk burning rate, 

v(t)= [ Tt ^y = v lf f( r) dxdy 



JD H K J D H 

which extends the notion of front speed. We derived lower bounds for long time averages of 
V(t) which behave like CA for large A, with constant C depending on the geometry of the flow. 
These bounds are valid for a class of flows that we call percolating. They are characterized 
by infinite tubes of streamlines connecting ±oo and include shear flows as a particular case. 
Our bounds imply the estimate ca > CA for traveling waves. Audoly, Berestycki and Pomeau 
gave a formal argument |]J suggesting that in the case when the shear flow varies on the scale 
much larger than the laminar flame width, one should have ca ~ A. One of the by-products 
of this paper is a rigorous proof of this conjecture. 

Our main goal in the present paper is to consider advection effects for a different physically 
interesting situation, where initial data are compactly supported. In this case, two generic 
scenarios are possible. If the support of the initial data is large enough, then two fronts form 
and propagate in opposite directions. Fluid advection speeds up the propagation, accelerating 
the burning. However, if the support of the initial data is small, then the advection exposes 
the initial hot region to diffusion which cools it below the ignition temperature, ultimately 
extinguishing the flame. 
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We consider for simplicity periodic boundary conditions 

(6) T(t,x,y)=T(t,x,y + H) 
in y and decay in x: 

(7) T(t,x,y) — > as |x| — > oo. 
We take 

H 

u(y)dy = 0. 



o 



A constant non-zero mean can be easily taken into account by translation. We consider 
the case when the width of the domain is larger than the laminar front width length scale: 
H 3> I = We will always assume that initial data T (x,y) is such that < T (x, y) < 1. 
Then we have < T < 1 for all t > and (x, y) G D. Moreover, we assume that for some L 
and 7] > we have 

(9) T (x,y) >9 + V for \x\ < L/2, 

T (x, y) = for \x\ > L. 

The main purpose of this paper is to study the possibility of quenching of flames by strong 
fluid advection in a model (Q). The phenomena associated with flame quenching are of great 
interest for physical, astrophysical and engineering applications. The problem of extinction 
and flame propagation in the mathematical model (|2|) under conditions (|3[), (0), was first 



studied by Kanel in one dimension and with no advection. He showed that, in the absence 
of fluid motion, there exist two length scales L < L\ such that the flame becomes extinct for 
L < Lq, and propagates for L > L±. More precisely, he has shown that there exist Lq and L\ 
such that 

(10) T(t, x, y) — > as t — > oo uniformly in D if L < Lq 

T(t, x, y) — > 1 as t — >• oo for all (x, y) e D if L > L 1 . 

In the absence of advection, the flame extinction is achieved by diffusion alone, given that the 
support of initial data is small compared to the scale of the laminar flame width I = k/vq. 
However, in many applications the quenching is the result of strong wind, intense fluid motion 
and operates on larger scales. There are few results available for such situations in the 
framework of the reaction-diffusion model. Kanel's result was extended to non-zero advection 
by shear flows by Roquejoffre who has shown that (|10|) holds also for m ^ with L Q and 



Li depending, in particular, on A and u(y). The second (propagation) part in ( |10"D was also 



proved in |21[ for general periodic flows. However the interesting question about more explicit 
quantitative dependence of L , L\ on A and u(y) remained open. Is it possible to quench the 
initial data that previously lead to an expanding solution by increasing A, but not changing 
the profile? How does this possibility depend on geometry of the profile u(y)7 Anyone who 
has tried to light a match in the wind has some intuition about this phenomenon. Yet, the 
mathematical answer turns out to be surprisingly subtle. In this paper, we also limit ourselves 
to shear flows. We are interested in the understanding the behavior of L and L\ for large 
A. The answer depends strongly on the geometry of the flow. In some cases the maximal 
extinction size grows L ~ A, and in others even the propagation size L\ remains finite as A 
goes to infinity. In the first case, we will say that u{y) is quenching. 
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Definition 1. We say that the profile u(y) is quenching if for any L and any initial data 
T (x, y) supported inside the interval [— L, L] x [0, H], there exists A such that solution of 
becomes extinct: 

T(t, x, y) — > as t — > oo uniformly in D 

for all A > Aq. We call the profile u(y) strongly quenching if the critical amplitude of advection 
Aq satisfies Aq < CL for some constant C(u, k,vq, H) (which has the dimension of inverse 
time). 

The key feature that distinguishes quenching from non-quenching velocities is the absence 
or presence of large enough flat parts in the profile u(y). 

Definition 2. We say that the profile u(y) G C°°[0,H] satisfies the H-condition if 
(11) there is no point y G [0,-ff], where all derivatives of u(y) vanish. 

The H-condition guarantees that the operator 

d d 2 
dx dy 2 

is hypoelliptic ||12|| . The study of existence of smooth fundamental solutions for such operators 
was initiated by Kolmogorov fl6f . Kolmogorov's work with u(y) = y served in part as a moti- 
vation for the fundamental result on characterization of hypoelliptic operators of Hormander 



12 1 . The hypoellipticity of the operator ( |12"D plays a key role in our considerations. Our first 



result is that the H-condition implies strong quenching. 

Theorem 1. Let f(T) be an ignition type nonlinearity , and let u G C°°[0, H] satisfy the 
H-condition. Then u(y) is strongly quenching. That means that there exists a constant 
C(u,n,v o,H) > that may depend on H, k, t> and u(y) but is independent of A such 
that 

(13) T(t,x,y) — > as t — > oo uniformly in D. 

This flame extinction occurs whenever the initial temperature To(x,y) is supported in an 
interval [—L,L] x [0, if], with L < C(u, k,Vo, H)A . 

The next result shows that a plateau on the order of the laminar front width 

(14) 1 = k/v 

in the profile u(y) prohibits quenching. (And therefore the conditions in Theorem [1] are 
natural.) 

Theorem 2. There exist universal constants Cq,C\ > 0, such that, if u(y) =u = const for 
y G [a — h, a + h] for some a G [0, H] and h > ho — CqI, then 

(15) T(t, x, y) — > 1 as t — >• oo 

uniformly on compact sets, for all A G M. This flame propagation occurs whenever the initial 
temperature Tq(x, y) satisfying (j^) is supported in an interval [-L, L] x [0, H] with L > Ciho. 

An interesting by-product of the proof of Theorem |2] is an estimate for the speed of traveling 
front solutions of @ when the shear flow varies slowly on the scale of the laminar flame width 
I. Let us define 



\u+\h = max < mm 

0<yo<H \ y£[y Q -ho,yo+h ] 
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with ho = CI given by Theorem and u + = max(«(?/), 0). 

Theorem 3. The speed of the traveling front c satisfies the upper and lower bounds 

A\u+\ho < c < A 1 1 w + 1 1 oo + v . 

The upper bound of Theorem |3| is contained in J7| (it is shown there for KPP type reaction, 
but this immediately implies the corresponding bound for ignition nonlinearity by a simple 
application of maximum principle). The left hand side is close to Hw+Hoo if u(y) is slowly 
varying on the scale ho- This agrees with the formal prediction of Audoly, Berestycki and 
Pomeau [0], and also (up to the addition of vo) with the results of Majda and Souganidis [JT7 
in the homogenization regime n — > 0. 

Unlike hypoellipticity, the quenching property is stable to small L°° perturbations: a small 
enough plateau (on the scale of the laminar front width /) does not stop quenching. 

Theorem 4. For every do > in (Qj there exists a constant B > such that, if a profile 
u(y) satisfies the H condition outside an interval y G [a — h, a + h] with h < hi — Bl, then it 
is strongly quenching. 

Moreover, the strongly quenching profiles are generic in the following sense: 

Theorem 5. The set of all strongly quenching shear flows u(y) contains a dense open set in 
C[0,H] (here C[0, H] is the space of continuous functions on [0,H]). 

In Section |] we show that all results on quenching, namely, Theorems [l], f| extend to the 
case of full system ([[]) with Le ^ 1. 

Finally, in the last section we prove that initial data of sufficiently small size (of the same 
order as in the case u(y) = 0) will be quenched by any shear flow Au(y). 

Acknowledgment. We thank N. Nadirashvili for fruitful discussions, and in particular for 
suggesting that we look for radially symmetric sub-solutions. We also thank G. Papanicolaou 
for stimulating discussions. This work partially supported by ASCI Flash Center at the 
University of Chicago. PC was supported by NSF grant DMS-9802611, AK was supported 
by NSF grant DMS-9801530, and LR was supported by NSF grant DMS-9971742. 

2. Quenching by a shear flow 
We prove Theorem [I] in this section. 

Proof of Theorem [7]. It suffices to show that there exists some time t Q such that 

(16) T(t , x, y) < 0o for all (x, y) G D. 

Then it follows from the maximum principle that T(t,x,y) < 8q for all t > to, and hence T 
satisfies the linear advection-diffusion equation 

T t + Au{y)T x = kAT. 

for t > t , which implies ([13]). We actually show that flTBD holds at t Q = k/vq for sufficiently 
large A. Recall that f(T) < T, and hence T(t,x,y) can be bounded from above using the 
maximum principle as follows: 

(17) T{t,x,y)<${t,x,y)eW*. 
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Here the function x,y) satisfies the linear problem 

(18) $t + Au{y)$ x = kA$ 

$(0,x,y) = T o (x,y) 
®(t,x,y) = <$>(t,x,y + H). 

Furthermore, we have 



/oo 
dzG{t,x- z)^{t,z,y) 
-oo 



with the function ^/(t, x,y) satisfying the degenerate parabolic equation 

(19) % + A U (y)ty X = K^yy 

^(0,x,y) = T (x,y) 
®(t,x,y) = V(t,x,y + H) 

and 

1 fx 2 
G(t,x) = exp 



We note that if u(y) satisfies the H-condition (|Tl|) then the diffusion process defined by ( |19|) 
has a unique smooth transition probability density. Indeed, the Lie algebra generated by the 
operators d y and d t + u{y)d x consists of vector fields of the form 

I"' I + ui ~ y) ^ u "( y) 7T> • • • ' u(n) ^7T> ■ ■ ■ 

oy ot ox ox ox ox 

which span R 2 if u(y) satisfies (|TTD. Then the theory of Hormander and the results 

of Ichihara and Kunita |]T3| imply that there exists a smooth transition probability density 
PA.(t,x,y,y') such that 

H 

V(t,x,y) = J dx' J dy'p A (t, x - x' , y, y') T (x', y'). 

R 



In particular, the function pAif) is uniformly bounded from above for any t > ||13|| . Then 
we have 

T(t,x,y) < e^l«m)\\L~ y < e^H^IU- < e^\\p A (t)\\ L ~ \\T oUly . 
It is straightforward to observe that 

PA{t, x,y,y) = —Po{t, —x, y, y ) 
A v 

with p being the transition probability density for ( |I9"D with A = t> . That is, po satisfies 

?El + v u ( ) dpo = r d2p ° 
dt dx dy 2 

pa(0,x,y,y') = 5(x)8(y - y') 
Po(t,x,y) =po(t,x,y + H). 

Therefore we obtain 

T(t,x,y) < e^'/^lboWH^JlToll^ < Ae^^\\p (t)\\ L ^LH 
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and in particular at time to = we have 



(20) T(t ,x,y) < C^|bo(*o)||x« LH < 9 

as long as 

A > c \\Voit Q )\\ LTy LH 

Theorem [l] follows from ( p0|) as explained in the beginning of this Section. □ 

We prove now Theorem f| that shows that a sufficiently small plateau in the profile u(y) is 
not an obstruction to quenching. 

Proof of Theorem (|. Let us define the set 

D r = lx [a — r,a + r] . 
As before, it suffices to show that solution of flLSD satisfies 



(21) $(t ,x,y)<^, t = ^ 



e 



and that is what we will do. First, we split the initial data for (pi]) into two parts: one 
supported on a strip D hl , containing the flat part of u(y), and another supported outside it. 
We will choose h\ = C\l > h such that any solution of ( |T9"D that is independent of x and 
with initial data supported inside will be small at time to — K / v o- The second part is 
supported away from the strip D^, where u(y) is flat. Therefore for a sufficiently small time 
it behaves like a solution of (|I^) with advection satisfying the H-condition. 
We choose hi as follows. Let 4>(t, y) be a periodic solution of 

<p t = (j) yy , 0(0,?/) = (f) (y), < (j) Q (y) < 1 

given by 

<j>(t,y) = J2M0)e 2lJ7Ty/H - 4Ktn2f/H2 . 

Then we have 

<i>(t,y)-- 

A simple estimate shows that 



H 



H 



4re 2 7r 2 7' 2 TJ n , 



and hence 
(22) 

Therefore we have 



(f>(to,y) 



H 



< C 



K 



<f){to,x,y) < ^ 
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as long as 

(23) 110(0, y)||£j <C t l, 

where C\ = 9qHI j '(10(Z + CH)). Let us choose hi so that (^) is automatically verified for 
initial data supported on D 2 i n : 

h < Cil/A 

with Ci as in (|23|). Let us assume that the width of the interval [a — h,a + h], on which -u(y) 
is constant satisfies 

- T" 

This condition determines the constant I? in the statement of Theorem |j. We may now split 
the initial data for ( |19|) as follows: 

2o(a;,2/) < Xo(y) + tpo(x,y)- 

Here the smooth function Xo(y) is supported in the interval [a — 2hi, a+2hi] while the function 
ipo(x, y) is supported outside the set [a — hi,a+hi\. Both of these functions satisfy in addition 
< Xo(y)ii } o{ x iy) — 1- Then the function <&(t,x,y) satisfies the inequality 

a?, 2/) < x(*, 2/) + V'(*> ^ 2/) 

with the functions x an d i> satisfying (|TjJ) with the initial data Xo an d respectively It 
follows from our choice of hi that 

x{to,y) < ^, 

so it remains only to estimate i/j(ta, x, y). We will do it separately for (x, y) inside and outside 
of the strip D hl / 2 . For the points (x,y) G D hl / 2 we have: 

H 

^{t,x,y) = J dx' J dy'p A (t,x - x' ,y,y')ij; (x' ,y') < J dx' J dy'p A (t,x - x' ,y,y') 

R R \y'-a\>hi 

<r{mt)\>^}<% 

for sufficiently small t. Here W(t) is the one- dimensional Brownian motion with diffusivity 
k, and V denotes probability with respect to it. Thus (|2T| ) holds inside D hl / 2 . 

In order to estimate the function ifj(t,x,y) outside D hl / 2 we introduce a profile u(y) that 
coincides with u(y) outside of the interval [a — (1 + S)h, a + (1 + S)h], <5 1, and satisfies the 
H-condition on the whole interval [0, if]. We define the process (X(t),Y(t)) by 

(24) dX(t) = u(Y(t))dt, dY(t) = \^dW(t), X(0) = x, Y(0) = y. 

Consider the stopping time r which is the first time when Y(t) enters the interval [a — (1 + 
5)h, a + (1 + 8)h], Then we have 

ip(t,x,y) < P x , y {(X(t),Y(t)) G supp Vo} = P x , y {(X(t),Y(t)) G supp ^ |r > t} P y (r > t) 
+P X:V {(X(t), Y(t)) G supp ^ |r < t} P y (r < t) 

< P X)lJ {{X{t),Y{t)) G supp ^ \r > t} Py(T >t) + Py(T < t). 
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Here P x<a denotes probability with respect to the process (X(t), Y(t)) starting at (x, y), while 
P y denotes probability with respect to Y(t) starting at y (recall that Y(t) is independent of 
x). The process X(t) for t < r is identical to the process (X(t),Y(t)) defined by (pl|) with 
u(Y) replaced by u(Y). Therefore we have 

if>(t,x,y) < P x , y {(X(t),Y(t)) E supp ^ |r > t} P y (r > t) + P y (r < t) 

< P x>y {(X(t), Y(i)) E supp ^0} + P y {r < t). 

Recall that (x, y) Dh x /2 and h < hi/ A. Therefore the point y is a fixed distance away from 
the interval [a — (1 + 8)h, a+ (1 + 8)h]. Hence we may choose t\ < to sufficiently small so that 



P y (r < h) < ^. 



10' 

Furthermore, the function ip(t,x,y) = P x<y <j (X(t),Y(t)) G supp ipQ [■ satisfies (|l^) with the 
initial data 



4>(0,x,y) 



1, (x, y) e supp tp 
0, (x, y) i supp 



However, u(y) is quenching and thus we may choose A so large that 





Therefore we have at t = t\: 



ip(t u x,y) < 



and hence the same upper bound holds at t = to > t%. Therefore ( pi]) holds also outside 

and Theorem || follows. The fact that u(y) is strongly quenching follows from this property 

of u{y). □ 

Theorem |5] is a simple corollary of Theorem |l|. 



Proof of Theorem |3|. The set of all profiles u(y) satisfying H-condition is dense in C[0, H], so 
by Theorem |I] the set of strongly quenching profiles is dense. To complete the proof, we will 
show that if u(y) satisfies H-condition, then there exists 5(u) such that if 

\\u(y) - u(y)\\ C [o,H] < S(u) 

then u(y) is strongly quenching. From the proof of Theorem |I|, we know that there exists 

a constant C(u) such that the solution ty(x,y,t) of the equation ( p!9|) with advection Au(y) 
satisfies 

*(x,y,t ) < - 



if the initial data ^f(x, y, 0) is supported on the interval [—L,L] with L < C{u)A (recall 
that such inequality implies quenching of T with the same initial data). Let ty(x,y,t) be a 



10 
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solution of flH3| ) with advection Au(y), and initial data supported in [— L', L']. Then by the 
Feynman-Kac formula (see, e.g. ||) 

t 

<f(x, y, t) < P x>y {x + Aju(y + W(s)) ds E [-L', L']) < 

o 

t 

P x>y (x + A j ' u(y + W(s)) ds e \-U - ASt, V + A8t\) < $(x, y, t), 



where fy(x,y,t) is the solution of (]TP|) with advection Au(y) and initial data equal to the 



characteristic function of the interval [—V — A5to, L' + A5to] for t < t . Now choose 5 < 
C{u)/t . Then 

6 

V(x,y,t ) < V(x,y,t ) < - 

e 

provided that L' < {C{u) — 5to)A, and hence u(y) is strongly quenching. □ 

3. Quenching for a system 

All results on quenching for equation (0) proved in the previous section extend directly 
to the case of the system ([j]). In this respect, the situation is similar to the case u(y) = 0, 



where all results on quenching proved by Kanel [14[] for a single equation extend to the system 
case. We make an assumption g(T) < T, which is just a normalization and corresponds to 
the condition f(T) < T in the case of a single equation. We take compactly supported initial 
data for the temperature, while for the concentration we assume < n(x,y,0) < 1. Notice 
that by the maximum principle, n(x, y, t) < 1 for every t. Therefore, T satisfies 

T t + u{y)T x - kAT < T, 

and from this point the analysis proceeds in the same way as for scalar equation (starting 
from (|TTD). We summarize the results in 

Theorem 6. Assume that g(T) = for T < 6q, and g{T) < T. Suppose that the velocity 
profile u(y) G C^°[0,H] satisfies H-condition, then it is strongly quenching. The same con- 
clusion holds if the profile u{y) has a plateau of the size h < h\ = Bl (where B is a universal 
constant). Moreover, the set of all quenching shear flows contains a dense open set in C[0, H}. 

4. Flame propagation 

We prove Theorem [| in this section. Let T(t, x, y) be solution of (0) with the initial data 
as in (0). We will use the following result of Xin [|2l| that holds for more general types of 



advection (its version for the shear flow was also proved by Roquejoffre in Consider 



(25) T t + u(x, y)-VT = kAT + —f(T) 

K 

T(0,x,y) = T (x,y)EL 2 (D) 
with u(x,y) being periodic in both variables, and < T (x,y) < 1. 
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Proposition 1 (Xin). Assume that the initial data in (fM) are such that 



lim T (x,y) = uniformly in D 

\x\^oo 

and 

T (x, y) > 9 + r] for \x\ < L. 
Then there exists Li(i],u) that depends on i] and u(x,y) such that if L > L\ then 

1, if ci < s < c r 



lim T(x — st, y, t) 



0, if ci > s or s > c r 



Here c\ and c r are the speeds of left and right traveling waves, respectively (ci < 0, c r > 0). 

The right traveling wave is described by ([|), (f|), and the left traveling wave satisfies (f|) 
and 

lim U(s,y) = 0, lim U{s,y) = 1, U.{s,y) > 0. 

s— »— oo s^+oo 

Notice that Theorems [I] and |] imply an estimate L\ > CA if u has no flat parts larger than 
certain critical size. On the other hand, Theorem § shows that L\ = h Q is independent of A 
if u has a sufficiently large flat part. 

Proof of Theorem f|. The proof of Theorem |2| proceeds in several steps. We consider the initial 
data satisfying ([5]). First we find h such that there exists a C 2 function (p(x,y) such that 
< < Oo + V, and 

v 2 

(26) «;A0+ > 

K 

and vanishes on the boundary of the disc of radius ho centered at the point (0, a) : 

0Uo=O, n = B((0,a);h ). 

Then in the system of coordinates that moves with the speed u the function <fi(x, y) is a sub- 
solution of (fj) in Qq- Therefore, initial data that start above <fi will not decay to zero. Next 
we consider the special solution $(t, x,y) of (||) with the initial data given by 

(27) *»(*,»)={*'*• " 

We show that x, y) satisfies (p~5| ) and that implies that ([15]) holds for arbitrary initial data 
2o ^ ^0) m particular, such as described in Theorem 0. 

Step 1. Construction of a stationary sub-solution. Choose 0±,02 so that 0q + i] > 
2 > 0i > O and define /i(T) by 

r 0, T <0i 

h{T) = \ m ^~ Bl \ e 1 <T<0 2 

{ f(T), 1 02 <T < I 

The function /(T) is Lipschitz continuous, and hence we may choose 0\ and 02 so that 
fi{T) < /(7 1 )- Therefore if satisfies 

(28) K A0+^A(0) = O 

K 
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then cj) satisfies (PB[). We are going to construct an explicit radial solution <p(r) of fl2"8| ) with 
the "initial" conditions 

0(0) = e 2 , ^(o) = o. 

Indeed, <p(r) is given explicitly by 

(29) 4>(r) = ft + (ft - ft) Jo ( ^) for r < ft - K6 



6*2 - 0i v oy /a 
Here Jo(0 is the Bessel function of order zero, and £i is its first zero. Furthermore, we have 

(30) <p(r) = Bin-, for R x <r 

R 

with B and R determined by matching and ( pCf ) at r = Then we get 

0i 



it = <4i W . exp 



/(*2 



CL I = -. 



Then d>(r) satisfies 



KA<p + —fU) > 0, 0(r) > 0, < r < R 

K 

<t>{R) = 0. 

Thus we will take the critical size of plateau in the velocity profile to be 2R so that the disc 
of radius R will fit in. 

Step 2. A sub-solution. Let us now assume that h > ho = R. We make a coordinate 
change 

£ = x — ut 

In new coordinates we have a function T(t,£,y) that solves 

(31) T t + A{-u + u(y))Ts = kA^T + ^/(T) 

with the initial data given by 

' <I>(Z,V), (^y)en = B((0,a);h 0) 



K 



Observe that 0(£, y) satisfies Q2"fi| ) inside Q since u(y) = u in fi . Moreover, T(t, £, y) > </>(£, ?/) 
on 9f2o, where vanishes. Therefore the maximum principle implies that inside Qq we have 

(32) T(t, C, y) > 4>(£, y) for all* > and (£, y) E Q . 
Note that Th(t, £, y) = T(t + h, £, y) solves (|3~ID with the initial data 

(33) T fc (0,e,!/) = T(A > e,l/)>r(0,e,J/). 

The inequality in ( j33"l) follows from (|3~2"D inside flo and the fact that T(t, x, y) > outside f2 . 
Therefore we have 

T{t + /i, f , y) - T{t, £, y) > for all /i > 0, t > and (f , y) E D 

and thus the limit 

T(£,y) = lim T(t,e,2/) 
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exists since T < 1. Moreover, the standard parabolic regularity implies that T(t,£,y) con- 
verges to T(£, y) uniformly on compact sets together with its derivatives up to the second 
order. Therefore T satisfies the stationary problem 

(34) (-u + u{y)% = kA^T + 
We also have 

T(£,2/)>0(£,2/)for (£,y)en . 
It is easy to show using the sliding method of Berestycki and Nirenberg || that for any 

(35) T{£,y) ></>{£ -r,y) 

where the right side is any translation of along the £ axis. Indeed, assume that there exists 
the smallest (say, positive) r such that T(£ ,2/o) = 0(£o — r , Ho) at some point (£o,yo) G Dh- 
Then the strong maximum principle implies that T(£, y) = 0(£ — r, y) for all (£, y) inside 
the translate f2 r = fi — re i of the disc Qq. But that contradicts the fact that T(£, y) > 
(/>(£ — r, y) = on the boundary dQ r . Then (^) implies that 

(36) f (£, a) > fl 2 = sup 0(f , y), for all (el. 

(£,v)en 

The next two lemmas show that ([34] ) and (|3~6"D imply that T(£, y) = 1. 
Lemma 1. Lei T be a solution of such that fiJBJ holds. Then we have 

(37) J d£dyf(f(x,y)) < oo, J «/|VT| 2 < oo. 

D D 

Proof. In order to show that integral of f(T) is finite we integrate (133) over the set (— L + 
C, L + C) x [0, H] with L large and ( e [0, /]. We get 

H H 

A y ds/[-w + u(j/)] [T (L + C,y)~ f(-L + (,y)] = K f dy%{L + {,y)- %{-L + C, y)} 
o o 

H L+C 

+^ jdy f de/(T(e,y) 

-L+C 

and average this equation in £ 6 [0, /]: 

/ H 

jjdcj dy[-u + u(y)][f(L + (,y)-f(-L + Cy)] 
o o 

H 

= jj dy[T(L + l,y)- T(L, y) - f(-L + l,y) + f(-L, y)] 
o 



+ v ljd(Jdy J de/(T(e,y). 



-L+C 
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Therefore we obtain 

H L 

iJ d V J dU{T{U))<— + 2HA[u+\\u\\ 00 } 

-L+l 

for all L, and hence the first inequality in fl37|) holds. In order to obtain the second inequality 
we multiply ([34]) by T and perform the same integration and averaging as above. This leads 
to 

l H 

dC I dy[-u + u(y)][f 2 (L + C,y)-f 2 (-L + C,y)} 



o o 

H 



K 

21 



dy[T\L + l,y)- f 2 (L, y) - T 2 {-L + l,y)+ f 2 (-L, y)] 



(38) 



I H L+C 

+± I d( I ' dy I di 



r„,2 



-T/(T) + «|VT| 



-L+C, 

and then the second inequality in ( j3?D follows from (|38|) and the first inequality in (|37|). □ 



Lemma 2. T/ie /zmzi function T(£, y) = 1 

Proof. Notice that T can not achieve local minima in D as follows from the maximum prin- 
ciple. Therefore if we define = min T(£, y) and fi(ot) = min T(a, y) then 

Q <£</3>0<y<-ff 0<y<H 

Ha, p = Ma) or /V/? = - Furthermore, if T(^,y) = 1 at some point, then T(£,y) = 1 
everywhere by the strong maximum principle. In particular, if /i(0) = 1, then T(^,y) = 1. 
Therefore we consider only the case that /z(0) < 1 and argue by contradiction. We have 
either /i(£) < /i(0) for any £ > 0, or for any £ < 0. Otherwise the minimum of T over the 
set [—£,£] x [0,H] would be achieved inside. Let us assume without loss of generality that 
/i(£) < /i(0) for any £ > 0. Consider 

. f i-fi(o) e 2 -e Q 

o = mm , 

V 2 ' 2 

For any £ > I, we have three options. 

1. T(£,y) G [yu(£),yu(£) + 5] for any y G [0, if], and > 6*o + 5. In this case, by definition 
of 5, properties @ of /, and since < /i(0) < 1, we have 



H 



f(T(£,y))dy>CH. 



2. < 6> + 5. Inequality ( |36] ) implies that there exist y%, y 2 such that T(£,yi) 
f(Z,y 2 ) = 6 2 . Then 



62-6, 



!/2 



f(T(£,y))dy 



> C\y 2 - 2/i | 
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2/2 / J/2 x 



J \T y \ 2 dy > \y 2 - yi\ x \j \T y \dy\ >5 2 \y 2 -yi\ 1 . 



2/1 \yi 



1)2 



Therefore, 

o 

(39) J f(f{t,y))dy J \f y \ 2 dy>C5 2 

H yi 

in this case. 

3. The remaining option is that > 9q + 5, and there exists y G [0, H] such that T(£, y) > 
MO + ^- I 11 this case, an argument identical to the reasoning of option two leads to the same 
bound (|39|). 

Overall, we see that for any £ < /, 



| y) + I Vf| 2 (£ )2/ )) dy>C, 



o 

where C depends only on H and 5. But this contradicts directly Lemma III. □ 

Using Proposition [I], we can now complete the proof of Theorem 0. Notice that even 
though we showed T(£, y) = 1, we still have to show the limiting function is the same in 
the original coordinates. Lemma implies that there exists a time t so that for all £ G 
[0, L (9 2 - O , Au(y))) and all y G [0, if] we have 

T{t ,Z,y) >e + v . 

Here L (rj, v) is defined by Proposition 0. Then we may apply Proposition |l| with s = in the 
original coordinates (x, y), and initial data T(t , x, y), and get (|15f) . The fact that q < < c r 
follows from @ and, for instance, results of 0. □ 

The sub-solution we constructed in Theorem is also useful for a proof of Theorem |||. As 
mentioned in the Introduction, we only need to show the lower bound, the upper bound is 
contained in [|7| . We begin the proof with two auxiliary lemmas. 

Lemma 3. Let u(y) be a function that is constant on an interval: 

u(y) — u, y G [a — h,a + h] 

that may or may not satisfy the mean zero condition (TJj. Assume that h > ho with ho given 
in Theorem^. Then we have ci(u) <u< c r (u) (recall that q(m) and c r (u) are the velocities 
of unique left and right traveling fronts). 

Proof. Let To be an initial data as in Theorem |2] such that 

o[ ,y) ~ \ o, (^y)^no 

with the function <p(x, y) constructed in the proof of Theorem ^| and given explicitly by ( p9[ ) 
and (pOl). Then we have 

T(t,x- ut,y) > <f>(x,y) 

and therefore Proposition |I] implies that q < u < c r because T(t, x — ut, y) may not go to 
zero as t — ► oo. □ 
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Let uq be a point at which |it+|/i is achieved: u(y) > \u+\h for y G [yo — h ,y + h ]. We 
define now a new velocity field 



(40) v(y) 
so that u(y) > v(y). 



u(y), y ^ [yo - h ,y + h ] 
|«+|ho, y G [yo - ^o,2/o + ho] 



Lemma 4. Lei u(y) > t>(y) 6e £u;o velocity fields. Then the corresponding right traveling 
front speeds satisfy c r {u) > c r {v). Similarly, if u(y) < v(y), then q(m) < Ci(v). 

Proof. Let T u,v (t,x,y) be solutions of 

v 2 

T, u + Au(y)T^ = kAT u + -^f{T u ) 

K 



and 



v 2 

T t v + Av(y)T v x = kAT v + ^f(T v ) 



K 

with the same initial data 

T u (0,x,y)=T v (0,x,y) = T (x,y). 

The function T (x, y) is assumed to be front- like and monotonic. That is, T (x, y) = 1 for 
x < —L, T (x, y) = for x > L and 

Then by maximum principle applied to the equation for dT u ' v /dx we have 

-(t,x,y) < for all t > 



dx 
and hence 

Tf + Au(y)T: - kAT v - V ^f{T") = A(u(y) - v(y))T" < 0. 

K 

Therefore by maximum principle we have 

(41) T u (t,x,y)>T v (t,x,y). 

However, by results of [21], |20| similar to Proposition [l] but for the front-like initial data we 
have 



and similarly for T u . Then c r (v ) > c r (u) would be incompatible with (41), and hence c r (v ) < 
c r (u). Naturally, an analogous result is true for left traveling fronts. □ 

Proof of Theorem [|. In order to finish the proof of Theorem [3] we apply Lemma [| to v (y) 
given by ([40]), and observe that c r {v) satisfies the estimate in Theorem |3| according to Lemma 

s □ 

Finally, we remark that in the situation of Theorem ^| strong advection not only does not 
quench the flame, but has directly opposite effect, according to Theorem |3| and Proposition [TJ. 
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Corollary 1. Let u(y) = u = const for y G [a — h, a + h] for some a G [0, H], and h is larger 
than or equal to the critical size h Q . Then, provided that the size of initial data L > h , we 
have 

1, if ci < s < c r 
0, if Q > s or s > c r 



lim T(x — st, y, t) 



Moreover, c r > \u + \h and c\ < — |w_|ft . 

Proof. The first statement is a direct corollary of Proposition [l| and Theorem ^. The second 
statement is the content of Theorem □ 

5. A LOWER BOUND FOR THE QUENCHING SIZE 

Recall that the burning rate is defined by 

dxdy 



v(t)= / T t (x, y y 



H 



D 



We also say that nonlinearity /(T) is of concave KPP class if 

/(0) = /(l) = 0, f(T) > for < T < 1, f"{T) < 0. 

We have previously shown [0 that for such nonlinearities and for front-like initial conditions 
the burning rate in the presence of any advection is bounded from below by Cv , more 
precisely: 

(42) V(t) >Cv (l-e- Cv ° t/K ). 

The physical meaning of is that no advection may slow up the burning significantly. We 
now show that similarly there is a fixed size of initial data that is quenched by any shear flow. 
That means that no shear flow may help prevent quenching of initial data with a fixed small 
support. Let T(t,x,y) be solution of 

(43) T t + u(y)T x = kAT + ^f(T) 

K 

with the initial data T(0,x,y) = T (x,y) and periodic boundary conditions @. 

Proposition 2. There exists a constant C > that depends only on the non-linearity f(T) 
such that given any initial data T with 



/T (x,y)dxdy < CI 2 , I = —, 



D 

we have 

(44) T(t, x, y) — > uniformly in D as t — > +oo 

for allu(y) G C[0,H}. 

Proof. We will show that solution of (|43l) with f(T) replaced by a KPP nonlinearity f(T) = 
MT(1 — T) and the same initial data drops below 6 before time to — K / v o- That will imply 
(H). Choose M such that f(T) < f(T) = MT(1 - T) (such M exists since / is Lipschitz 
continuous) and let T be the solution of 

f t + u{y)f x = KAf +^/(T) 

K 
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with initial data To and periodic boundary conditions (|6|). Let us also define 

V(t) = [ f t (x,y) d -^ = M [ f(l-f) dxdy 



H k J ' H 

b d 



Note that T(t, x, y) > T(t, x, y) and hence Proposition [| follows from the following Lemma. 

Lemma 5. There exists a constant C > that depends only on the non-linearity f(T) such 
that given any initial data Tq with 



/T (x,y)dxdy < CI 2 , I = — , 
v 



D 

and any u(x, y) G C 1 (D) , there exists a time t\ < to such that T(t,x,y) < 6q for all (x, y) 6 D. 
Proof. The following inequality holds for V(t): 

V(t) + ^^>2n [\Vf\^. 
V ; Mvl dt ~ J 1 1 H 

D 

We define the set 

S(t) = {ye [0,H] : 3x e R such that T(t,x,y) > 9 /e 2M } C [0,H). 

Observe that there exists a constant 5 > that depends only on M and 8q such that if the 
Lebesgue measure |«S(r)| < 51 at some time r e [0, to/2], then ||T(to)|U°° — ^o- Indeed, we 
have then 

f{r,x,y)< X (y)+ip{x,y), \\x(v)\\<Sl, ip{x,y) < 6 /e 2M . 
Therefore we have with C as in (|22|) 



T(t ,x,y) < e 



M 



(C + 1)S 



e 2AI 



< e n 



for a sufficiently small 5. Hence it suffices to consider the case when \S(t)\ > 51 for all 
t G [0, to/2]. We claim that then at any time t G [0,to/2] we have 

(45) J \VT\ 2 dxdy J f{f)dxdy > CI 2 

D D 

with the constant C depending only on M and do- This may be proved similarly to Lemma 
2 in [[/J. Therefore we obtain 

~ _^dV{t)_ > 2Cv 2 l 2 
U Mv 2 dt ~ H 2 V(t)' 

Thus we have 



Therefore 



V(t) > for t /2 <t<t . 
H 



T{to,x,y)> J T (x,y)dxdy + Cnt . 



D D 
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However, we have an a priori bound 

f(t , X ,y)d X dy<e-jT„ {x ,y)d X dy 

D D 

which contradicts the previous inequality if 

(46) J T (x,y)dxdy < J^l^ = C M l 2 - 

D 

Therefore T(t,x,y) has to drop below 9 if the initial data satisfies ( [46] ) and Lemma || is 
proved. □ 

Proposition [| then follows from Lemma |5|. □ 

Remark. The uniform quenching size in Proposition |2| is optimal. Indeed, the subsolution 
we have constructed in Section [| for shear flows with a flat part has L l norm C7 2 , and thus 
one cannot expect that initial data with L 1 -norm larger than CI 2 are quenched by all shear 
flows. 
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